Orientational orders in binary mixtures of hard HGO molecules 
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Based on a standard constant-NPT Monte Carlo molecular simulation, we have studied liquid crystal 
phases of binary mixtures of non-spherical molecules. The components of the mixtures are two kinds 
of hard Gaussian overlap (HGO) molecules, one kind of molecules with a small molecular-elongation 
parameter (small HGO molecules) cannot form stable liquid crystal phase in bulk, and other with 
a large elongation parameter (large HGO molecules) can form liquid crystal phase easily. In the 
mixtures, like the large HGO molecules, the small HGO molecules can also form an orientation- 
ordered phase, which is because that the large HGO molecules can form complex confining surfaces 
to induce the alignment of the small molecules and generate an isotropic-anisotropic phase transition 
in the whole binary mixtures. We also study the transition on different mixtures composed of 
small and large HGO molecules with different elongations and different concentrations of the large 
molecules. The obtained result implies that small anisotropic molecules might show liquid crystal 
behavior in confinement. 



I. INTRODUCTION 

Surface-induced ordering of liquid crystal (LC) 
molecules is of important technological and scientific in- 
terests [1,2]. Very recently, Boamfa [3] et al. experimen- 
tally observed an isotropic-nematic (IN) surface phase 
transition in a mixture of nematic LCs on a substrate. 
Usually, in confined LCs, interaction between the LC 
molecules and the surfaces of the confining walls gives rise 
to a nematic phase at the surface, though the bulk is still 
an isotropic phase [3-5]. Based on Monte-Carlo (MC) 
simulation methods, Gruhn and Schoen [6,7] have studied 
the microscopic structure of molecularly thin confined LC 
films. They found that the orientationally ordered thin 
LC films were in thermodynamically equilibrium with 
the isotropic bulk phase. In our recent work [8], based 
on MC simulations of ellipsoid-like molecules confined in 
slit pores, we found a similar phenomenon in confined 
LC molecules. Additionally, as ellipsoid molecules with 
a small molecular-elongation parameter are confined in 
very thin slit pores, the surface of the pores can induce 
an orientation-ordered phase which cannot be formed in 
bulk due to the very small elongation parameter of the 
studied molecules. The result implies that the surface of 
confining walls not only induces surface IN phase transi- 
tion in LC materials, but also stabilizes the orientation- 
ordered phases in non-LC materials (in the bulk, non-LC 
materials such as small-molecular-weight fluids cannot 
form aniostropic phases). Since small- weight molecules 
can response very fast to the variation of external fields, 
the small-molecular-weight LCs would be of great advan- 
tage in many applications. However, in simple confined 
geometry, the surface-induced effects are limited to the 
several molecular layers neighboring the surfaces of the 
confining walls [3,8], but in complex confined geometry, 
we may expect surface-induced effects in whole fluids. 

Polymer-dispersed liquids can be used to study the 
properties of fluids confined in complex geometry, since 



the dispersed polymers form complex networks. Re- 
cently, Chiccoli et al. studied the network-induced or- 
dering in LCs based on MC simulations [9]. In the mix- 
ture of polymers and LC molecules, the polymers are 
thought as "surfaces" to confine the LC molecules and 
the surface effect exists in the whole system. We may 
address the following question: Do the found properties 
near the confining walls in simple confined geometry glob- 
ally exist in the complex confined geometry (i.-e. mix- 
tures)? However, we are interested in a particular ques- 
tion in this paper: Since small ellipsoid-like molecules 
(no orientation-order phases in bulk) can be induced an 
IN transition in a thin slit pore [8], can we find simi- 
lar results in complex pores? For example, in mixtures 
of the small ellipsoid-like molecules and other molecules, 
can IN transition occur? We will consider a binary mix- 
ture composed of small ellipsoid molecules and usual LC 
molecules. These LC molecules can form complex sur- 
faces to affect the small molecules. Can the surfaces in- 
duce ordering phases in the small (non-LC) molecules? 
Since there are large amount of small-weight molecules 
in the mixture, the systems may show some characteris- 
tics of expected small-molecule LCs. 

In this paper, based on standard constant-pressure 
(NPT-ensemble) Monte Carlo (MC) simulations, we 
study binary mixtures of hard Gaussian overlap (HGO) 
molecules. One component of the mixtures is large HGO 
molecules with large molecular elongation parameter k 
(defined as the length-to-breadth ratio), and other is 
small HGO molecules with small k. The large HGO 
molecules mimic LC molecules which can form LC phases 
in bulk alone, whereas the small HGO moleucles, corre- 
sponding to normal small anisotropic molecules, cannot 
form stable LC phases in bulk. Hence each mixture is a 
solution with LC solute and non-LC solvent. Mixtures 
of LC molecules and other molecules have been studied, 
such as earlier researches of LCs with small amounts of 
impurities [10-12] and recent studies on mixture of poly- 
mers and LCs [9,13]. But few works focused on the or- 
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dered phases of the solution of LCs in small molecules. 
Here, we calculate the equation of states and orienta- 
tional ordering parameter of the mixtures. Similar to our 
previous finding that confining walls induce an alignment 
of small HGO molecules in thin slit pores [8], we find, in 
the mixture, the large-fc (LC) molecules can induce an 
alignment of the small-A: (non-LC) molecules, hence form 
an orientation-ordered phase in the whole system. Our 
results indicate that small molecules will form LCs with 
the help of LC molecules in their mixtures. 

The paper is organized as follows: In Sec. II, we 
describe a HGO model and give an explicit formula 
of intermolccular interaction between different-fc HGO 
molecules, then we describe the presentation of our simu- 
lations. Simulation results are shown in Sec. III. Finally 
we conclude and discuss in Sec. IV. 



II. POTENTIAL MODEL AND 
COMPUTATIONAL DETAILS 

There is a great number of molecular simulations in the 
literature on the properties of LCs. Among these sim- 
ulations, LC molecules are described using anisotropic 
interactive models [7,14-23]. Although the detail re- 
sults depend on the models used, some general conclu- 
sions have been found: Repulsive interactions between 
molecules play key roles to the LC structures, so some 
hard nonspherical models [14-20] have been widely em- 
ployed to study the properties of LCs; LC phases can 
form only when the anisotropic parameters of molecular 
interactions are greater than the critical values. Other- 
wise (i.e. for molecules with small anisotropic interac- 
tive parameters), as one increases the pressure or de- 
creases the temperature, these systems will freeze or 
crystallize before forming LC phases [8,24]. Many re- 
cent works focus on the macroscopic properties and the 
phase transitions of models with large anisotropic param- 
eters [6,7,18,19,21,25]. 

Among LC models, the HGO model [8,14-16,24] is 
widely used. The model came from the Gaussian overlap 
potentials developed by Berne and Pcchukas [26]. It is 
supposed that the mass (or electron) density of molecules 
was a Gaussian function of space vector x, 



Gj(x) = exp(-x • A i 1 • x), 
A = (erf? - u^UiUi + all. 



(1) 



Here Uj is a unit vector along the principal axis of the ith 
molecule and I is the unit matrix, <7ii and crj_ correspond 
to the size of the molecules along and perpendicular to 
its principal axis, respectively. The interaction between 
two molecules i and j is supposed to be proportional to 
the overlap of the two molecules, 



V(ui,Uj,rij) ~ |Aj 



-1/2, 



"7 



dxG l (x)G J (x - Tij) 



A, + A, |- 1/2 exp[-r • (A, + A^)" 1 • r] 



e(uj, Uf) cxp[-r 2 /cr 2 (u J , Uj, f)], 



(2) 



where, r and f are the length and the unit vector of the 
vector rij, respectively. We have, 



cr 2 (u i ,u J ,f) = f ■ (Ai + Aj) 1 -f. 



(3) 



If the three eigenvalues and corresponding unit eigenvec- 
tors of Ai + Aj are noted as Xk and Vk, {k = 1,2,3), 
respectively, we have 



a 2 (u 4 ,Uj,f) = ^(f • v k f/\k- 



(4) 



fc=i 



Actually, since a vector perpendicular to and Uj is an 
eigenvector of Aj + Aj, we only need to diagonalize a 
2x2 matrix. Since ik (k — 1,2,3) are unit vectors and 
perpendicular each other, we have X^fc=i(T ' ^k) 2 = L 
Using these relationships, we easily obtain, 

v(u u Uj ,r) = a {l - f [f^(r • + f^(r • v 2 ) 2 ]}^ 2 

and 6 = v /«?. + (l-4-)("i-u J ) 2 . 
en as, 

'-r, -<*l) + (4-«3l) 



where 
c = 



Xij and aij are writen as 



Xij — 



(°?II+<*l) + (°ji 
(^l-^x)-(^l 



°5±V 



+ K„-^)' 
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respectively. The directions of the unit eigenvectors v\ 
and V2 are given as, 



Vl 

V2 



)uj + (uj • Uj)(l - aij)uj, 



{Ui ■ Uj)(l 



Ui 



(7) 
(8) 



respectively. Obviously, if the ith molecule and the jth 
molecule are identical, we obtain the usual formula of 
identical molecules [26], 

X (f ■ Uj + f ■ Uj) 2 (f ■ Ui — f ■ UjY 



a(Ui,Uj,r,) = a {l - — [ % — ^— + % — 

2 1 + X (Ui ■ Uj) 1 — X(Ui-Uj) 

where X^ is simply noted as X. 

The interactive potential in HGO model is simplified 

as, 
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V(Ui,Uj,Tij) 



if 



< (j(Ui,Uj,fi 



0, if Tij > a(ui,Uj,fij). 



(10) 



where a(u i: Uj,fij) is defined as Eq. (5). The HGO 
model is thought as an approximation for the excluded 
volume of hard ellipsoids of revolution [16,24]. In this 
paper, we study binary mixtures consisting of two kinds 
of HGO molecules with molecular elongation kh and ks, 
respectively. Here the elongation k is defined as, 
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h = <J l \\/cn±- (11) 

Here i = L or S, corresond to large and small HGO 
molecules, respectively. For simplification, wc suppose 
(Tl± = as±, then a is a constant which is set as unit of 
length in our simulations, and can be writen as 

Xij = (kj - 1 + k] - \)/{kj + 1 + k] + 1). (12) 

We use standard constant-pressure (const-NPT) MC 
simulations to obtain the orientational orders and the 
equation of state P(p) ( where P corresponds to the real 
Pco/ksT). The simulations are performed with N = 500 
or 256 molecules which consist of Nl large-elongation 
molecules and N — Nl small-elongation molecules. The 
selected ks is about 2.0 and k^ is between 3.0 and 
6.0. According to a large number of previous simula- 
tions [8,14-16,24], the selected large-A: molecules can form 
orientation-ordered phases in bulk alone, but the selected 
small-fc molecules cannot form stable LC phases in bulk 
up to a typical high density of liquids [8,24]. In our 
simulations of the mixtures, the simulating box is cu- 
bic and equally fluctuates in three directions, cubic peri- 
odic boundary conditions are used. The simulations are 
organized in MC cycles, each MC cycle consisting (on 
average) of N trial translational and rotational molecu- 
lar displacements and one trial volume fluctuation. The 
trial displacements of two kinds of molecules are sepa- 
rately treated. The maximum step length of each trial 
move is automatically chosen at each pressure for mak- 
ing the acceptable probability fall between 0.4 and 0.5. 
The starting configuration is a face-center-cubic lattice 
which melt at low pressure and equilibrate for 10 5 MC 
cycles. The system is slowly compressed in small pres- 
sure steps. For any given pressure, the system is typically 
equilibrated for 10 5 MC cycles and an average is taken 
over all configurations of additional 10 5 MC cycles, then 
the final configuration is set as the starting configuration 
of the next pressure. At some higher pressure (density) 
zones, 3 ~ 5-times longer MC cycles are used in equili- 
brated and averaged processes. 

The orientational order parameter S is calculated in 
the simulations as the largest eigenvalue of the ordering 
Q a p tensor, defined in terms of the components of the 
unit vector u ia along the principal axis of the molecules, 

1 N 3 1 
Qa/3 = ^2(-u ia u il3 - -8 a[j ) (13) 

i=l 

Besides detecting the order parameter of all molecules, 
we also separately calculate the order parameter (Sl 
and Ss) of larger and smaller molecules by limiting the 
average of all molecules in Eq. (13) to single kind of 
molecules. 



III. RESULTS 

We first simulate mixtures composed of Nl large 
HGO molecules (k L = 5) and N s = N - N L small 



HGO molecules (ks = 2.0). The concentration of large 
molecules is x — Nl/N. We present the equation of 
state and orientation order parameter 5* of the mix- 
tures with different x in Fig. 1. In label of the figures, 
N L /N is listed. For pure small HGO molecules(fc = 2.0), 
the simulations are performed with N = 256 molecules, 
but for other cases, we simulate systems with N = 500 
molecules. In Fig. 1 (a), the equations of state exhibit a 
discontinuity for all larger-x systems. From Fig. 1(b), we 
can clearly find the spontaneous phase transition between 
isotropic and anisotropic phases in the HGO system con- 
posed of small amounts of large-fc molecules. However, 
the order parameter S of the mixture is smaller than 
unity, which implies that the molecules only partially 
align or only partial molecules align. With decreasing 
concentration of large molecules x, the transition pres- 
sure increases and the corresponding order parameter S 
in the anisotropic phase decreases. 

To detect the microscopic structure of two kinds of 
molecules in the mixtures, we calculate separately the 
orientation order parameter of large and small HGO 
molecules. In Fig. 2, we show relations between the 
three order parameters S and the pressure in a mixture, 
where x = 0.4, ks = 2.2 and fc^ = 5.0. These order 
parameters describe the alignment of large-fc molecules, 
small-fc molecules, and all molecules, respectively. From 
Fig. 2, we find a transition for both large and small 
molecules beginning at a same pressure P ~ 2.0 in the 
mixture. With increasing pressure, the large molecules 
almost completely align, but the small molecules only 
partial align. Thus the average S of all molecules is far 
smaller than unity. We also find the direction of the large 
molecules is parallel to that of the small molecules (not 
shown). Therefore, it can be concluded that the align- 
ment of the small molecules is unambiguously ascribed to 
the alignment of of the large molecules. Additionally, the 
transition seems to occur at a pressure range, rather than 
at an exact value, which may correspond to the known 
I-N coexistence in LC mixtures. In Fig. 3, we show the 
orientation orders of large molecules and small molecules 
with different fci and ks values. The IN transition value 
of small molecules nematic order parameter decreases as 
the elongation of large molecules fci increases. This sur- 
prising result is ascribed to the lower transition pressure 
as kL is larger (see Fig. 3(d) ). 

IV. CONCLUDING REMARKS 

From a great number of researches of LCs, people had 
known that only molecules with larger anisotropic pa- 
rameters can form LC phases in bulk. For example, the 
molecular elongation of normal LC molecules is larger 
than 4. In our previous work [8], we found that small 
anisotropic molecules may form stable LC phases as the 
molecules are confined in very thin pores. The result 
indicates that some normal inorganic small molecules 
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(whose k is about 2) may form LC phases in confined 
geometries (or near surfaces). In current widely stud- 
ied polymer-dispersed LC systems, the polymers can be 
thought as complex surfaces confining LC molecules, so 
we guess that small anisotropic molecules may also show 
LC phases due to the existing large (LC or polymer) 
molecules in the mixtures. 

In this paper, we present results from standard NPT- 
ensemble simulations of mixtures of small and large HGO 
molecules. Among the two components, the small HGO 
molecules cannot form LC phases alone in bulk, but 
large HGO molecules can. We detect possible isotropic- 
anisotropic phase transition in the mixtures. The results 
presented here give clear evidence of spontaneous for- 
mation of orientation-order phases in the mixtures. We 
find that the spontaneous ordering of the large and small 
molecules occur at the same pressure. It implies that the 
ordering of the small molecules is ascribed to the induce- 
ment of the large molecules. We also study the induced 
phase transition behaviour for different largc-molccule 
concentrations and different molecular elongations of the 
large and small molecules. 

The confinement-induced LC phases in mixtures 
present here may be useful. In the mixtures, the viscos- 
ity should be smaller than large-molecule systems (usual 
LCs), so we think that both the large and the small 
molecules may faster align in external field, which implies 
that the mixtures may show the characteristics of ex- 
pected small-molecule LCs, such as short response time. 
However, it is necessary to study the confined effects 
based on more reliable molecular models. 
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Fig.1 X. Zhou et al. k s = 2.0 and k L = 5.0, 10 A 5 MC 
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FIG. 1. (a) The Equation of state (pressure versus den- 
sity) of binary mixtures with two kinds of HGO molecules. 
Their molecular elongation k are 5.0 and 2.0, respectively. 
The data in label correspond to Nl/N of each curves. For 
Nl = 0, we select N = 256. (b) The Orientation order pa- 
rameter S versus pressure of the mixtures. All parameters 
and symbols are same as that in Fig. 1. It is clearly shown 
an Isotropic-nematic phase transition except in small x cases. 
However, except a trivial case (a; = 500/500 = 1.0), in nc- 
matic phase, the order parameter 5* of the mixture does not 
arrive at unit. 



Fig.2 X. Zhou et al,k s =2.2 and k L =5.0, x=200/500 
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FIG. 2. Orientation order parameters of large, small and 
all molecules versus pressure P in a mixtures. Where 
x = 200/500, k s = 2.2 and k L = 5.0. 



FIG. 3. Orientation order parameters of large and small 
molecules versus pressure in mixtures with different ks/kL 
for detecting the elongational dependence of the ordering pa- 
rameters. For systems with fci = 3.0, Nl/N = 101/256; for 
other systems, N L /N = 200/500. The symbols of (b) are 
same as that of (a) and the symbols of (c) are same as that of 
(d). (a) Orientational order parameter Sl of large molecules 
in mixtures versus P (only near transition zone); (b) Orien- 
tational order parameter Ss of small molecules in mixtures 
versus P; (c) Sl versus P; (d) Ss versus P. In (c) and (d), 
for the case of ks/kL = 2.0/5.0, we only show data in lower 
pressure range (triangle). 
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Fig. 3 X. Zhou etal. x=0.4 
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